PASJ: Publ. Astron. Soc. Japan , 1-??, 
(c) 2008. Astronomical Society of Japan. 



Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: II. The 

Molecular Gas Disk 

Hiroyuki NAKANISHI 1 and Yoshiaki SOFUE 2 
l Nobeyama Radio Observatory, Minamimaki, Minamisaku, Nagano 384-1305 
2 lnstitute of Astronomy, The University of Tokyo, 2-21-1 Osawa, Mitaka,Tokyo 181-0015 

hnakanis @ nro. nao. ac.jp 

(Received 2006 February 14; accepted 2006 August 3) 
Abstract 

We created a three-dimensional distribution map of molecular gas throughout the Milky Way Galaxy using the 
latest 12 CO( J =1-0) survey data cube and rotation curve based on the kinematic distance. The radial distribution 
of the molecular gas shows a central peak and a second peak around 0.5 i?o- The thickness of the molecular disk 
slightly increases from 48 pc to 160 pc with the galactocentric distance within a radius range of 0-1 1 kpc. We were 
able to trace the Outer, the Perseus, the Sagittarius-Carina, the Scutum-Crux, and the Norma arms as logarithmic 
spiral arms with pitch angles of 11° — 15°. Considering that the pitch angles of the spiral arms are within this 
range, the Norma and the Outer arms seem to be identified as the same spiral arm. We could also trace a midplane 
displacement, whose amplitude is nearly constant inside a 10 kpc radius and increases beyond this radius. The ridges 
of midplane displacement form leading spiral arms. 
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1. Introduction 

The global three-dimensional (3D) structure of the Milky 
Way Galaxy is one of the essential informations to study the 
Galactic objects in detail. In Paper I (Nakanishi & Sofue 2003), 
we have analyzed the latest 21 -cm line survey data, and pre- 
sented the 3D distribution map of H I gas in the Milky Way 
Galaxy. In this paper, we deal with the molecular gas distribu- 
tion as Part II of this series. 

A survey in the 12 CO(J =1-0) line is a good tool for trac- 
ing the interstellar matter (ISM) in the inner disk of the Milky 
Way Galaxy (inside the solar orbit), while the H I data is ad- 
vantageous for study of ISM in the outer disk (outside the so- 
lar orbit). One of the most famous CO line surveys is the 
Columbia CO survey, whose earlier survey was carried out with 
the 1 .2-m telescopes in New York City, USA, and Cerro Tololo, 
Chile (Dame et al. 1987; Bronfman et al. 1988; Bronfman et 
al. 1989). This survey covered the entire Milky Way with a 
grid spacing of 0?5. Another higher resolution survey is the 
Massachusetts-Stony Brook Galactic plane CO survey, which 
was made using the 14-m telescope of the Five College Radio 
Astronomy Observatory (FCRAO) (Sanders et al. 1986). The 
covered ranges of the galactic longitude and latitude were from 
I = 8° to 90° and from b = -l.°05 to +1P05 in the first quad- 
rant, respectively. The grid spacing was 3' x 3' or 6' x 6'. These 
large CO surveys were ideal for exploring into the molecular 
gas distribution in the Milky Way Galaxy. There have been 
several studies on the Galactic structure using these data. 

Dame et al. (1986) and Myers et al. (1986) showed face- 
on views of the first Galactic quadrant examining positions 
of the largest individual molecular complexes identified in the 
Columbia CO survey of the first Galactic quadrant (Cohen, 
Dame, & Thaddeus 1986). They found clear evidence for the 



Sagittarius spiral arm. Cohen et al. (1985), Grabelsky et al. 

(1987) , and Grabelsky et al. (1988) showed distribution of gi- 
ant molecular clouds in the fourth quadrant and they clearly 
delineated the Carina arm. Clemens, Sanders, & Scoville 

(1988) presented the molecular gas distribution in the first 
quadrant of the inner disk of the Milky Way Galaxy using 
the Massachusetts-Stoneybrook CO survey, found a ring-like 
structure at the radius of 0.5i?o (Ro- the Galactocentric dis- 
tance of the Sun), and traced the Sagittarius and Perseus arms. 
The distance of the molecular gas was determined by the kine- 
matic distance using a rotation curve. The near-far problem 
in the inner Galaxy was overcome by introducing a vertical 
model of the molecular gas distribution. Solomon & Rivolo 

(1989) showed a face-on view of the Milky Way Galaxy ex- 
amining individual giant molecular clouds identified in the 
Massachusetts-Stoneybrook CO survey. 

Hunter et al. (1997) presented the molecular gas distribu- 
tions in the entire Galactic disk of the Milky Way Galaxy using 
the Columbia survey (Dame et al. 1987) for modeling the dif- 
fuse gamma-ray emission. The distance of the molecular gas 
was determined by the kinematic distance. Although there is 
a near-far problem in the inner disk, they divided gas density 
equally at near and far points for 6 = 0°. 

Sawada et al. (2004) determined the molecular gas distribu- 
tion in the Galactic center by considering the absorption of the 
OH line to determine the line-of-sight distance of the molecu- 
lar gas, using the CO data cube of Bitran et al. (1997). They 
found that molecular gas is elongated in the line-of-sight direc- 
tion against the Sun, and that this feature is associated with a 
central galactic bar. 

Thus, there are several studies on the molecular gas distribu- 
tion. However, there has been no three-dimensional gas den- 
sity map of the whole Galaxy. Dame, Hartmann, & Thaddeus 
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(2001) published a new compilation of Columbia 12 CO( J = 1- 
0) survey of the whole of the Milky Way Galaxy. The sampling 
and the sensitivity of this survey were much improved over that 
of the old Columbia survey. This is a great tool for investigat- 
ing the molecular gas distribution throughout the entire Milky 
Way Galaxy. In this paper, we show a molecular gas distribu- 
tion in the Milky Way Galaxy using the Columbia 12 CO( J =1- 
0) survey data and the latest rotation curve. 

2. DATA 

We used the latest compilation of 12 CO(J =1-0) survey 
data (Dame, Hartmann, & Thaddeus 2001), which were ob- 
tained with the 1.2-m telescopes in New York City, USA, and 
Cerro Tololo, Chile. The half-power beam width (HPBW) was 
8.' 4, which corresponded to 2.4 pc at the heliocentric distance 
of 1 kpc. The data consisted of 12 CO( J = 1-0) line spectra over 
the entire Milky Way. We used data of surveys named 2, 8, 18, 
31, 33, and 36. Grid spacing of these data was 0P125 - 0?25. 
The total radial velocity (V T ) coverage was 332 km s _1 and ve- 
locity resolution was 1.3 km s -1 . The r.m.s. noise level per 
channel was 0.12 - 0.43 K. We used spectra within a Galactic 
latitude range of \b\ < 1.5°. 

The rotation curve, which we used to transform the radial 
velocity of the CO gas into the distance, is the same as we used 
in Paper I. The inner and outer rotation curves are taken from 
Clemens (1985) and Dehnen and Binney (1998), respectively. 
The Galactic constants i?o = 8.0 kpc (the Galactocentric dis- 
tance of the Sun) and Vq = 217 km s _1 (the solar rotational 
velocity) are adopted. 

3. METHOD 

3.1. Binding the Data 

3.1.1. Inner Galaxy 

For the inner Galaxy, we averaged the CO data over five 
channels (> 6.5 km s _1 velocity width) for each spectrum 
in order to eliminate the cloud-cloud motion of 3 - 3.9 km 
s _1 (Clemens 1985; Alvarez et al. 1990) and to obtain a higher 
signal-to-noise ratio (S/N ratio). 

The molecular gas is distributed in a more clumpy form than 
the H I gas. Since we intend to develop our understanding of 
the global structure of the Milky Way Galaxy, we must smooth 
the CO emission on a larger scale in order to eliminate such 
a clumpy features and to model the vertical structure of the 
molecular gas disk. Since the typical scale of the random mo- 
tion of the molecular clouds is 0.015i?o kpc (Clemens 1985), 
we averaged the data over a 240 pc width when the line-of-sight 
length corresponding to 6.5 km s _1 velocity width was smaller 
than 240 pc. Thus the CO data were averaged over > 6.5 km 
s _1 velocity width and > 240 pc width. 

In addition, we averaged the data in the longitudinal direc- 
tion to eliminate the clumpy features of CO clouds and to ob- 
tain a much higher S/N ratio. For the first (0° < / < 90°) and 
fourth (270° < I < 360°) quadrants, we bound the data, so that 
the longitudinal width at the tangential point is equal to 240 
pc-width. 

Thus the CO data were averaged over > 240 pc in the line-of- 
sight direction and longitudinal direction except for the inner 
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region of the tangential radius, where the data were averaged 
over > 240 pc width only in the line-of-sight direction. The 
data were not averaged in the latitudinal direction at all. 
3.1.2. Outer Galaxy 

For the outer Galaxy, we averaged the data in the longitu- 
dinal direction so that longitudinal separation becomes 0?5, 
which is nearly the same as the H I data. Moreover, we av- 
eraged the CO data over two channels (> 2.6 km s _1 velocity 
width), which is also the same as the H I data. 

Figure 1 shows the points where the molecular gas densities 
were calculated. 

3.2. Molecular Gas Distribution in the Outer Galaxy 

In this paper, we choose a Cartesian coordinate whose x— 
axis coincides with the line crossing the Sun and the Galactic 
center and whose origin coincides with the Galactic center. The 
z-axis is chosen so that it is parallel to the rotational axis. The 
Sun is located at (—8 kpc, 0,0). We also choose a cylindrical 
coordinate (R,0,z) so that the angle 6 = 180° coincides with 
the direction toward the Sun and the angle 6 = 90° is parallel 
toZ = 90°. 

The distance of the molecular gas was determined by kine- 
matic distance in this paper, the same way as in Paper I. The 
distance of molecular gas in the outer disk (R>Rq) is uniquely 
determined, and then we obtained the molecular gas distribu- 
tion in the outer disk using the equation 

AV r 

nn 2 =XT h —— (1) 
Ar 

where nn 2 denotes the volume number density of the molec- 
ular hydrogen H2, X the CO-to-P^ conversion factor, and r 
the heliocentric distance. We adopted two kinds of conver- 
sion factor; (1) a constant conversion factor of X — 1.8 x 10 20 
(H 2 cm" 2 K" 1 km" 1 s) (Dame, Hartmann, & Thaddeus 2001) 
and (2) a varying conversion factor with the Galactocentric dis- 
tance (Arimoto, Sofue, & Tsujimoto 1996). Arimoto, Sofue, & 
Tsujimoto (1996) calculated X using the data in Solomon et al. 
(1987). 

Since the conversion factor can be approximately expressed 
by exponential function of the Galactocentric radius, we fitted 
an exponential function to the resultant radial profile and ob- 
tained the following equation, 

X[H 2 cm- 2 K- 1 km- 1 s] = 1.4exp (iJ/llkpc). (2) 

3.3. Molecular Gas Distribution in the Inner Galaxy 

There is, on the other hand, a near-far ambiguity in the inner 
disk (R > R ); i.e. there are two heliocentric distances satisfy- 
ing the same radial velocity. Therefore, we could not determine 
the H2 distribution in the inner Galaxy in the same way as in 
the outer Galaxy. In order to divide molecular distributions at 
near and far points, we introduced the vertical structure of the 
Galactic disk. Following Spitzer (1942), we assumed that the 
vertical distribution of the molecular gas is expressed with 

n(z) = n sech 2 (£), (3) 

e = log(l + V2)i^, (4) 

Zl/2 
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where z and 21/2 denote the height from the Galactic plane 
(b = 0°) and the scale-height, respectively. At z = z 1 / 2 the 
density, n, gets a half value of the maximum at the midplane. 
The subscript means values at the midplane. The scale-height 
was determined at the tangential points in the inner Galaxy, 
which is described in the next subsection. 

Assuming that the vertical distribution of the H2 gas follows 
equation (3), the H2 column density, 7Vh 2 (6), is represented 
as a function of the Galactic latitude b using four parameters: 
™h 2 i> ri H 2 2^ z 0l , an d zq 2 . That is, 



iV H2 (&) = n H2 isech 2 (£;i 



An 

cos b 



where 



a = iog(i 



£ 2 = log(l + V2) 



r\ tan b — 

Zl/2 

Ti tan b — zq 2 

Zl/2 



n H2 2sech 2 (^ 2 )^4, (5) 
coso 



(6) 



(7) 



For each point shown in figure 1, we sought a parameter set 
consisting of 71h 2 i> ^h 2 2> zq ± , and zq 2 , which gives the small- 
est residual in subtracting the model distribution from the ob- 
servation. The fitting resolution of the four parameters is 0.01 
H2 cm~ 3 for njj 2 i and Uh 2 2> an d 1 pc for Zq 1 and z$ 2 , respec- 
tively. We restricted zq 1 and zq 2 within a range of —250 pc to 
+250 pc. The grid spacing is 0P125, which corresponds < 17 
pc at tangential points, at most points in the inner Galaxy. Since 
there are enough points (more than three points) to resolve near 
vertical distribution, it is possible to divide the total CO emis- 
sion into near and far emissions. We thus obtained parameter 
sets of 71h 2 i, hh 2 2, zq x , and z$ 2 for discrete points in the inner 
Galaxy. 

To show how well the fitting method works, we present 
a residual map in figure 2. In order to make this map, we 
first measured standard deviation of residuals between a fit- 
ted model curve and observed data. And then, we calculated 
relative residuals, which is a fraction of the standard deviation 
relative to the peak of the observational data. This residual map 
shows that the relative residual is less than 20 - 30 % in most 
of the region, and that the fitting method does not work well 
near the Sun. 

3.4. Scale-height 

Tangential points which satisfy R = Rq sinZ present no near- 
far problem in the inner Galaxy. Hence, the vertical structure 
at the tangential points can be determined uniquely. 

The gas of velocity range |V^ max — cr | < \V r \ < 00 is lo- 
cated at the tangential points, where V T , V r ,max, an d °" are ra- 
dial velocity, maximum radial velocity, and random velocity, 
respectively. The random velocity was taken to be 5 km s _1 
(Clemens 1985). We thus obtained vertical distributions, and 
fitted a model distribution p(z) = posech 2 (£) to determine a 
scale-height of z 1 / 2 - Figure 3 shows a full-width half max- 
imum (FWHM) (2z 1 / 2 ) plotted versus the Galactocentric ra- 
dius, which is a mean value in a 1 kpc bin. For a comparison, 
figure 3 shows also FWHM of molecular disk calculated by 
Bronfman et al. (1988) and Grabelsky et al. (1987), who use 
the same data that we use. We adopt the Galactic constant Rq 
of 8 kpc although Bronfman et al. (1988) and Grabelsky et al. 



(1987) use R of 10 kpc. Therefore, the radius and FWHM 
are reduced by factor of 8/10. Filled squares, open circles, and 
open triangles denote the data from this study, Bronfman et al. 

(1988) and Grabelsky et al. (1987), respectively. Our result is 
almost consistent with Bronfman et al. (1988) and Grabelsky 
etal. (1987). 

Figure 4 shows comparisons of a model and observations. 
A vertical distributions of molecular gas can be approximately 
reproduced by equation (3). 

— Figure 3 — 

3. 5. Construction of Three-Dimensional Map 

We thus obtained parameter sets of the molecular distribu- 
tion at each point of figure 1 . 

First, the parameter sets obtained for the discrete points 
in the inner Galaxy were averaged with a Gaussian function 
whose FWHM is 0.48 kpc, in order to eliminate data with bad 
fitting results. Then we calculated the H2 distribution in the 
inner Galaxy using the averaged parameter sets. The H2 gas 
densities obtained at the discrete points in the outer Galaxy 
were transformed into continuous distribution by interpolat- 
ing with the same Gaussian function. Averaging length for the 
z-direction was fixed at 35 pc, which corresponds to the grid 
spacing of 0P25 at the heliocentric distance of 16 kpc. Finally, 
we combined the H2 distributions in the inner and outer Galaxy 
to obtain a CO distribution cube. In this paper, the grid spacing 
of the cube is fixed at 200 pc for the x and y directions and at 
20 pc for the z directions, and the size of the resultant cube is 
40 kpc x 40 kpc x 2 kpc in x, y, and z direction, respectively. 

Since the number of points where the molecular gas densities 
were calculated (figure 1) in the azimuthal range \9\ < 30° is 
too few to make a map, we do not show this region. 

— Figure 1 — 

4. Results 

4.1. A Face-on Map 

Figure 5 shows a face-on map of the molecular H2 gas in the 
Milky Way Galaxy. 

Figure 6 shows a face-on map of the molecular gas when the 
conversion factor is calculated with equation (2). In this map, 
the molecular gas distribution in the outer region is slightly 
enhanced. However, there is no large difference between these 
maps. Hence, our discussions below about the structure of the 
molecular gas distribution are based on the H2 map obtained 
adopting a constant conversion factor for simplicity. 

— Figure 5 — 



— Figure 6 - 



Figure 14 ■ 
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4.2. Vertical Sliced Maps 

Figures 7 and 8 present vertical cross sections of the resultant 
cube sliced with planes which are parallel to the z-axis. We 
made 6 slices at the azimuthal angles in each figure. Panels of 
azimuth angles of 6 = 40°, 60°, 80°, 100°, 120°, and 140° in 
figure 7 are in exactly the opposite directions: 9 = 220°, 240°, 
260°, 280°, 300°, and 320° in figure 8. 

— Figure 7 — 

— Figure 8 — 

5. Structure of the Molecular Disk 

5.1. Radial Distribution 

Figure 9 shows a radial distribution of the surface density av- 
eraged azimuthally over the 1 kpc-width radius when the con- 
stant conversion factor was adopted. For a comparison, fig- 
ure 9 also shows radial distributions calculated by Bronfman 
et al. (1988) and Grabelsky et al. (1987). The radius is re- 
duced by factor of 8/10 because of difference in the adopted 
Galactic constant. Moreover, we adopted the conversion factor 
of 1.8 x 10 20 although Bronfman et al. (1988) and Grabelsky 
et al. (1987) used X of 2.8 x 10 20 . Hence, we reduced the sur- 
face density by factor of 1.8/2.8. The data from Bronfman et 
al. (1988) is a mean value of the northern and southern data. 
Filled squares, open circles, and open triangles denote the data 
from this study, Bronfman et al. (1988) and Grabelsky et al. 
(1987), respectively. The surface density decreases from 26.7 
Mq pc -2 to 0.3 Mq pc -2 with the Galactocentric distance. 
Our result is almost consistent with Bronfman et al. (1988) and 
Grabelsky et al. (1987). The plotted values are given in table 1. 

Radial distributions of the molecular gas in spiral galaxies 
are often expressed with exponential functions (Young et al. 
1995). However, this figure shows that the radial profile is 
not expressed by an exponential function. The surface density 
is quite high at the Galactic center and near the 4 kpc radius, 
which is in response to the central condensation and Four-kpc 
ring found in the face-on map. The Galaxy is known to have 
a bar in the central region (e.g., Nakada et al. 1991). These 
double peaks are thought to relate to the barred structure (e.g., 
Nakai 1992). 

— Table 1 — 

— Figure 9 — 

We also show the radial distribution of the molecular gas 
density at the midplane in figure 10. Radial distributions de- 
rived by Bronfman et al. (1988) and Grabelsky et al. (1987) are 
plotted for a comparison. The radius is reduced by factor of 
8/10, and the density is multiplied by (1.8/2.8) x (10/8) 2 be- 
cause of the difference in the adopted Galactic constant and the 
conversion factor. Filled squares, open circles, and open tri- 
angles denote the data from this study, Bronfman et al. (1988) 
and Grabelsky et al. (1987), respectively. The radial distribu- 
tion is similar to that of the surface density. The density de- 
creases from 8.73 cm -3 to 0.03 cm~ 3 with the Galactocentric 



distance. Our result is almost consistent with Bronfman et al. 
(1988) and Grabelsky et al. (1987). The radial profile of this 
value is not expressed by an exponential function, either. The 
plotted values are given in table 1 . 

— Figure 10 — 

5.2. Scale-height 

Figure 3 shows that the FWHM of the molecular gas disk 
slightly increases from 48 pc to 160 pc with the Galactocentric 
radius in the radius range of 0-11 kpc. The scale-height is 
smaller than that of the H I disk. The plotted values are given 
in table 1 . 

5.3. Midplane Displacement 

Figures 7 and 8 show that the midplane of the inner disk is 
slightly displaced from the Galactic plane, which is known as 
a midplane displacement, a tilted disk, or warped disk. The 
warping of the outer disk was clearly delineated by May et al. 
(1997). 

Figure 1 1 shows the midplane displacement versus the az- 
imuthal angle 9 relative to the Galactic center for individual 
radii. The vertical axis denotes the distance between the mid- 
plane of the molecular disk and the plane of b = 0°. We calcu- 
lated the midplane displacement, z m , by averaging z with the 
weight of density, n(z), 

J zn(z)dz 

Zrn — 7- 7 ■ (8) 

J n(z)az 

The midplane displacement can be expressed with a sinu- 
soidal function. Hence, the molecular gas disk can be con- 
sidered as a collection of tilted rings. We fitted the following 
sinusoidal function to the obtained data, 

z rn = z a + Asm{9-9 ). (9) 

The fitting algorithm used in this paper is the nonlinear least- 
squares Marquardt-Levenberg algorithm used in a software 
gnuplot. The obtained parameters are given in table 2. Figure 
12 shows that amplitude of the midplane displacement is nearly 
constant inside a radius of 10 kpc. However, the amplitude 
steeply increases beyond the radius of 10 kpc, where the H I 
warp is remarkable (Nakanishi & Sofue 2003). 

Figure 13 shows loci of the maximum points of individual 
tilted rings. The solid and dashed lines denote the loci where 
the midplane is displaced upward and downward, respectively. 
The loci form leading spiral arms, which is consistent with the 
rules of the warp summarized by Briggs (1990), and which is 
interpreted as follows: The midplane displacement can be con- 
sidered as a wave, called a bending wave (Binney & Tremaine 
1987). Defining fl and v as the rotational and vertical frequen- 
cies, respectively, the pattern speed of the bending wave is ex- 
pressed as Q, + v/m or O — v/m (to: integer). Since the first 
one denotes a fast propagating wave and a trailing arm, it dis- 
appears rapidly as a result of the winding process. On the other 
hand, the leading spiral component is long-lived. 

— Table 2 — 
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5.4. Central Concentration 

A highly concentrated molecular gas, which is called the 
central molecular zone (CMZ; Morris & Serabyn 1996), is seen 
at the Galactic center. The gaseous column density reaches 
> 25M Q pc -2 at the center when we adopt a conversion factor 
of 1.8 x 10 20 cm -2 K _1 km s _1 . Though it is plausible that 
abundant molecular gas is distributed near the Galactic center, 
the shape shown in this map is suspicious, because the non- 
circular motion is noticeable in the central region due to the 
central Galactic bar. Sawada et al. (2004) created a face-on 
view of the CMZ using the data from the CO (J =1-0) sur- 
vey of the Galactic central region Bitran et al. (1997). Because 
they derived it using without using the kinematic distance, their 
map is more reliable. The molecular gas in the CMZ is elon- 
gated in the line-of-sight direction. This feature is thought to 
be associated with a central galactic bar. 

5.5. Spiral Arms 

5.5.1. Scutum-Crux Arm 

Positions of the edges of the Scutum-Crux arm in the first 
and fourth quadrants are I = 33° and I = 305° (Georgelin 
and Georgelin 1976). In the direction of I = 305°, an arm- 
like structure can be detected around (R, 6) ~ (7kpc, 250°), 
which corresponds to a part of the Scutum-Crux. In figure 14, 
a schematic tracer of the Scutum-Crux arm is presented by a 
logarithmic spiral arm starting at (R,6) = (3.0 kpc,30°) with a 
pitch angle of 12°. The parameters of these logarithmic spiral 
arms used in figure 14 are given in table 3. 

— Table 3 — 

5.5.2. Sagittarius-Carina Arm 

The Sagittarius-Carina arm is a famous arm located near the 
Sun and is known to arc in the first and fourth quadrants (e.g., 
Georgelin & Georgelin, 1976; Cohen et al. 1985; Dame et al. 
1986; Grabelsky et al. 1987; Grabelsky et al. 1988; Clemens, 
1988; Solomon & Rivolo, 1989). According to Georgelin and 
Georgelin (1976), apparent edges of this arm are I = 50° and 
I = 283°. This arm can be traced by a logarithmic spiral arm 
with a pitch angle of 11°. In figure 14, a schematic tracer of 
the Sagittarius-Carina arm is presented by a logarithmic spiral 
arm starting at (R,6) = (2.0 kpc, 180°) with a pitch angle of 
11°. 

5.5.3. Perseus Arm 

A clear spiral structure can be found to arc from (R, 9) ~ 
(7kpc,90°) to (R,6) ~ (11 kpc,170°). This spiral arm is the 
so-called Perseus arm. This arm goes by the Sun outside the 
Solar orbit. The shape of this arm in our map is consistent with 
the result of Clemens, Sanders, & Scoville (1988). A cross sec- 
tion of this arm is found at R = 10 — 11 kpc in the top panel 
of figure 7. This arm can be traced more clearly in I — v dia- 
gram of Dame, Hartmann, & Thaddeus (2001) and is detected 
from I ~ 55° in the negative velocity to I ~ 270° in the positive 
velocity. 

This arm traced in the molecular gas map is consistent with 
that detected in the H I map (Paper I). It is more prominent 
in the H I map. The pitch angle of this arm was estimated to 
be 15°. In figure 14, a schematic tracer of the Perseus arm 
is presented by a logarithmic spiral arm starting at (R, 9) = 
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(5.7 kpc, 30°) with a pitch angle of 15°. 

5.5.4. Norma- Outer Arm 

At larger radius than the Perseus arm, the Outer arm can be 
traced from (R, 9) ~ (9 kpc, 50°) to (R, 6) ~ (11 kpc, 120°). 
The Outer arm can be seen also in the longitude-velocity dia- 
gram of the CO survey (Dame, Hartmann, & Thaddeus 2001). 
This arm appears in the Galactic longitude range of 30 < I < 
90°. The Outer arm is also detected in the H I map similarly to 
the Perseus arm. This arm traced in the molecular gas map is 
also consistent with that detected in the H I map, and is more 
prominent in the H I map. This arm can be traced with a loga- 
rithmic spiral arm with the pitch angle of 15°. 

The Outer arm can be found in the vertical sliced map (figure 
8). The Outer arm is located at the radius of R ~ 9 — 12 kpc 
in the 9 = 80°, 120°, and 140° panels of figure 8. Particularly 
in the 9 = 120° and 140° panels, this arm is displaced from the 
Galactic plane due to warping. 

The Norma arm is known to exist in the fourth quadrant 
(e.g., Georgelin & Georgelin, 1976). Georgelin and Georgelin 
(1976) stated that the apparent edge of this arm is / = 327°. 
In the fourth quadrant of our H2 gas map, there is an elon- 
gated lump with a high H2 column density in the direction 
of I = 327°, which corresponds to the Norma arm. Though 
Georgelin and Georgelin (1976) mentioned that the Norma arm 
is symmetrical to the Sagittarius-Carina arm, the Norma arm 
seems symmetrical to the Scutum-Crux arm rather than the 
Sagittarius-Carina arm. 

These two spiral arms seem to be the same arm, assum- 
ing the pitch angles of these arms are 15°. Hence, we call 
these arms the Norma-Outer arm. There are other possibil- 
ities of connections of (i) the Norma and Perseus arms and 
(ii) the Scutum-Crux and Outer arms. However, the radius of 
the Perseus is expected to be too small to be connected to the 
Norma arm in the fourth quadrant. Similarly, the radius of the 
Scutum-Crux arm is too large to be connected to the Outer arm 
in the first quadrant. Therefore, the connection of the Norma 
and Outer arms seems most natural. 

In figure 14, a schematic tracer of the Norma-Outer arm 
is presented by a logarithmic spiral arm starting at (R, 9) = 
(4.0 kpc, 240°) with a pitch angle of 15°. 

5.5.5. Comparison with arms identified from HII region data 
In order to evaluate the validity of the spiral arm that we 

trace, we compare our result with Georgelin and Georgelin 
(1976). We overlay the spiral arms that Georgelin and 
Georgelin (1976) traced on our tracing spiral arms in figure 
15. The sizes of spiral arms are changed so that the distance 
between the Galactic center and the Sun is the same. The 
Sagittarius-Carina, the Scutum-Crux, and the Norma arms that 
we trace are coincident with those traced by Georgelin and 
Georgelin (1976), although the radius of the Perseus arm is 
slightly different for I > 90° in the cases of our result and 
Georgelin and Georgelin (1976). 

5.6. Ring-like structure 

In the first quadrant (0° < I < 90°) a ring-like feature is found 
at i?~ 4 kpc, which corresponds to the ring-like structure found 
by Clemens, Sanders, & Scoville (1988). The difference in 
the scale of this structure is due to the adopted Rq. Clemens, 
Sanders, & Scoville (1988) adopted Rq = 10 kpc, while we 



3-D Distribution of Molecular Gas in the Milky Way Galaxy 



6 



H. Nakanishi and Y. Sofue 



[Vol. , 



adopted Ro = 8 kpc. Hence, the scale of our map is 4/5 times 
as large as the map of Clemens, Sanders, & Scoville (1988). 
We call this structure the Four-kpc Ring instead of the Five- 
kpc ring in this paper. 

The edge of the Four-kpc ring corresponds to I = 30° — 40°. 
This coincides with the edge of the Scutum-Crux arm. In ad- 
dition, the Sagittarius-Carina arm seems to be merged into the 
Four-kpc Ring in the azimuthal range of 9 < 120°. The Four- 
kpc ring may therefore be a complex of two spiral arms and 
remarkably enhanced as if it was a ring. 

Bronfman et al. (1988) pointed out that the radial distribution 
of molecular gas is flatter in the south hemisphere than in the 
north, and that the density peak at 0.5i?o in the north broadens 
in the south into a plateau. This statement is consistent with 
our conclusion that two spiral arms in the south are merged in 
the north. 

6. Possible Errors in our Molecular Gas Map 

Since the molecular gas map we obtained in this study is 
based on several assumptions, we must note the systematic er- 
rors due to deviations from these assumptions. 

First, we assumed that the Milky Way galaxy rotates ab- 
solutely circularly around the Galactic center. However, the 
gaseous motion deviates from circular motion due to the den- 
sity wave or the central Galactic bar, believed to exist at the 
Galactic center. The isovelocity contours are typically dis- 
torted in 1 kpc according to observations of external galaxies. 
Uncertainties concerning the location of the gas are about 1 
kpc. 

Second, we used the model expressed by equation (3) to di- 
vide the mixture emission from near and far points. Though 
the resultant map can change due to the adopted function, the 
global feature is almost the same for any function as shown in 
Paper I. 

Third, fitting method does not necessarily work. Typical 
residual is less than 20 - 30 % of midplane density (figure 2). 
However, residual becomes much larger near the Sun, which 
implies that the fitting method does not work well in this re- 
gion. 

Fourth, conversion factor also affects the resultant gaseous 
map. We tried to make 3-D maps using a constant conversion 
factor and a varying conversion factor with the Galactocentric 
distance. Outer structure is enhanced when we use the varying 
conversion factor. 

7. SUMMARY 

We constructed a three-dimensional molecular gas map of 
the Milky Way Galaxy using the new compilation data set 
of 12 CO(J =1-0) survey obtained with 1.2-m telescopes and 
the latest rotation curve, based on the kinematic-distance. We 
made molecular gas maps by adopting two types of conver- 
sion factor; one is constant, and the other varies with the 
Galactocentric distance. Global features are the same for both 
though the contrast is slightly different. 

We present a new 3-D map of the molecular gas in this pa- 
per. The resultant maps confirm the following features of the 
Galactic molecular disk. 



1 . The radial distribution of the molecular gas surface den- 
sity has double peaks at the center and at the radius of 
0.5Ro, which is a feature similar to barred spiral galax- 
ies. 

2. The FWHM of the molecular disk slightly increases 
from 48 pc to 160 pc with the Galactocentric distance 
within the radius range of 0-1 1 kpc. 

3. There is a midplane displacement, whose amplitude is 
nearly constant inside the radius of 10 kpc, increasing 
beyond this radius. The ridges of midplane displacement 
form leading spiral arms. 

4. Molecular gas arms corresponding to the Outer, the 
Perseus, the Sagittarius-Carina, the Scutum-Crux, and 
the Norma arms could be traced as logarithmic spiral 
arms with the pitch angle of 11° — 15°. Supposing that 
the pitch angles of the spiral arms are within this range, 
the Norma can be identified as being the same as the 
Outer arm. 

We are grateful to Dr. Thomas Dame et al. for providing us 
with their large CO survey data. 
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Fig. 1. Grid points where gas densities were calculated. The hashed area denotes a region where there are too few points to make a density map. 



Table 1. Radial variations of parameters of the molecular disk 



Radius 


Surface Density 


Midplane Density 


FWHM 


(kpc) 


(M Q pc- 2 ) 


(cm -3 ) 


(PC) 


0-1 


26.7±20.3 


8.73±6.91 


48±20 


1-2 


6.0±5.4 


1.28±1.11 


56±18 


2-3 


3.5±1.8 


0.76±0.42 


68±14 


3-4 


3.5±1.6 


0.75±0.32 


88±24 


4-5 


4.6±2.0 


0.83±0.38 


78±12 


5-6 


3.9±2.0 


0.66±0.38 


102±10 


6-7 


2.8±1.4 


0.48±0.27 


84±22 


7-8 


1.9±1.0 


0.26±0.17 


90±30 


8-9 


0.9±0.7 


0.09±0.07 


186±20 


9-10 


0.5±0.4 


0.05±0.05 


182±20 


10-11 


0.3±0.2 


0.03±0.03 


160±20 
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Fig. 2. Residual map. We here define the 'relative residual' to be a fraction of standard deviation of residuals between fitted model curves and observation 
data relative to the peaks of the observation data. 
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Fig. 3. FWHM of the H2 gas layer versus the Galactocentric radius. Filled squares, open circles, and open triangles denote the data from this study, 
Bronfman et al. (1988) and Grabelsky et al. (1987), respectively. 



Table 2. Fitted parameters of the midplane 



Radius 


zo 


A 


Oo 


(kpc) 


(PC) 


(PC) 


(°) 


0-1 


5±1 


13±1 


292±3 


1-2 


11±1 


36±2 


288±4 


2-3 


-2±4 


22±4 


265±15 


3-4 


-8±5 


12±6 


244±32 


4-5 


2±3 


23±4 


277±13 


5-6 


-12±3 


18±5 


190±13 


6-7 


-6±6 


10±9 


193±43 


7-8 


-9±6 


21±9 


133±23 


8-9 


-27±7 


40±9 


109±15 


9-10 


-49±10 


83±12 


76±11 


10-11 


-26±30 


160±25 


85±11 
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Fig. 4. Vertical distributions of the molecular gas. Horizontal and vertical axes denote height from the Galactic plane (6 = 0°) and H2 density, respectively. 
Thick and thin lines indicate observation data and model curves, respectively. 



Table 3. Parameters of schematic spiral arms 



Name 


Pitch Angle 


Beginning Point 


Endin, 


I Point 




(°) 


(kpc) 


(°) 


(kpc) 


(°) 


Scutum-Crux 


12 


3.0 


30 


9.1 


330 


Sagittarius-Carina 


11 


2.0 


180 


10.5 


310 


Perseus 


15 


5.7 


30 


12.6 


200 


Norma-Outer 


15 


4.0 


240 


12.3 


120 
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Fig. 5. Face-on molecular hydrogen map adopting the constant conversion factor of X = 1.8 X 10 20 (Dame, Hartmann, & Thaddeus 2001). The contour 
levels are 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, and 12.8 M Q pc~ 2 . 
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Fig. 6. Face-on molecular hydrogen map adopting the conversion factor X taken from Arimoto, Sofue, & Tsujimoto (1996), where X increases ex- 
ponentially with the galactocentric distance. The contour levels are the same as those of figure 5; 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, and 12.8 Mq pc~ 2 . 
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Fig. 7. Vertical cross-section maps. Each map shows the H2 volume densities in a sheet through the Galactic center perpendicular to the Galactic plane. 
First plate: Individual panel shows sections through 9 = 220°, 8 = 240°, 6 = 260°, 6 = 280°, 9 = 300°, and 9 = 320°. The contours are drawn at 
levels of 0.05, 0.1, 0.2, 0.4, 0.8, 1.6, and 3.2 H 2 cm- 3 . 




Fig. 8. Vertical cross-section maps. Each map shows the H2 volume densities in a sheet through the Galactic center perpendicular to the Galactic plane. 
Individual panel shows sections through 9 = 40°, 9 = 60°, 6 = 80°. = 100°, 6 = 120°, and = 140°. The contours are drawn at levels of 0.05, 0.1, 
0.2, 0.4, 0.8, 1.6, and 3.2 H 2 cm" 3 . 
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Fig. 9. Radial profiles of the H2 gas surface density. The vertical axis is a logarithmic scale. Filled squares, open circles, and open triangles denote the 
data from this study, Bronfman et al. (1988) and Grabelsky et al. (1987), respectively. 
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Fig. 10. Radial profiles of the H2 gas density at the midplane. The vertical axis is a logarithmic scale. Filled squares, open circles, and open triangles 
denote the data from this study, Bronfman et al. (1988) and Grabelsky et al. (1987), respectively. 
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Fig. 11. Midplane displacements for individual radii. Horizontal and Vertical axes denote the azimuthal angle around the Galactic center and the distance 
from the plane of b = 0° . 
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Fig. 12. The amplitude of the midplane displacement versus the Galactocentric distance. 
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Fig. 13. Loci where the midplane displacements reach the maximum in the face-on view. The thick and dashed lines denote the loci where the midplane 
is displaced upward and downward, respectively. 
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Fig. 14. Schematic tracers of H2 gas arms superimposed on the obtained H2 column density map; No.l: Norma-Outer arm; No. 2: Perseus arm; No. 3: 
Sagittarius-Carina arm; No. 4; Scutum-Crux arm. The contour levels are the same as those of figure 5; 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, and 12.8 Mq pc — 2 . 
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Fig. 15. Comparison of our result and spiral structure presented by Georgelin and Georgelin (1976). Thicker lines denote spiral arms which we traced. 
Thinner dashed lines denote spiral arms which Georgelin and Georgelin (1976) traced. 



